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Recently equations for the generating functional in the perturbative quantum chro- 
modynamics (QCD) have been extended by including the non-perturbative dissi- 
pation in QCD jets. The resulting equations have been solved rigorously and new 
family of scaling solutions , the so-called S - scaling, generalizing the well-known 
Kubo-Nielsen-Olesen scaling law for hadron multiplicity distributions have been 
found. The relevance of 8- scaling is discussed in the Landau - Ginzburg theory 
of phase transitions. Preliminary application of these ideas to the pp data of the 
UA5 Collaboration is presented. 



1 Introduction 

A suitable framework for the comparison of multiplicity distributions at dif- 
ferent energies is provided by the idea of KNO scaling EJ, which states that 
at sufficiently high energies : lim s _ >00 (< n > On)Q^ n a n)~ l = f( z ), where 
<7„ is the cross-section for producing exactly n particles, < n > is the average 
multiplicity of produced particles and f(z) is an energy-independent scaling 
function in z = n/< n >. The KNQ scaling is fullfilled in e + e~ collisions up 
to the highest energies consideredel . The e + e~ data are-well accounted for 
by the Perturbative Quantum Chromodynamics (PQCD) a calculations aad, 
in particular, by the Fragmentation - Inactivation - Binary (FIB) process oB, 
which is equivalent to PQCD with dissipation effects in jet cascading! . In pp 
collisions, KNO scaling holds up to ISR energies but, at higher .energies, UA5 
Collaboration reported significant breaking of the KNO scaling!! . 

The KNO scaling was derived from the hypothesis of Feynman scaling 
It is now known that Feynman scaling is violated at very high energies 
so other arguments should be looked for. In this context, a possible relation 
between the KNO scaling and the phase transitioaJn Feynman- Wilson gasQ as 
well as the criticality of self-similar FIB process crB was pointed out. It seems 
therefore that KNO scaling may have more profound reasons to appear. Since 
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new data on multiplicity distributions at still higher energies will be available 
soon , we believe that the discussion of information contained in multiplicity 
distributions is particularly urgent and challenging. 

2 Order parameter fluctuations 

Let us consider a self-similar statistical system of size N, in which phase 
changes are characterized by the order parameter 77. Self-similarity in such 
sytem means that for different scales which are characterized by different val- 
ues of order parameter < 77 >, fluctuations of normalized order parameter 
77/ < rj > are identical. We know that this is a feature of thermodynamical 
systems such as the Ising model at the critical point of second order phase 
transition. Defining the anomalous dimension as : 

g = lim JV^oo^^ (In < N\r)\ >) , (1) 

one can show that the partition function : Zn ~ N^ 1 - 1 ^ 9 ^ ~< \rj\ > and the 
probability density P[rj\ obeys the scaling law : 

< M > pm = = $ (f^y , (2) 

in the centered variable Zm = (77 — rj*)/ < \rj\ >, with 77* the most probable 
value of rj for given parameters of the system. We call (|^) the first scaling law. 
The scaling limit is defined by the asymptotic behaviour of P[rj\ when rj — > oo, 
< 1 77 1 >^ 00, but Z(i) has a finite value. The first scaling law. holds also for 
fluctuations of any power 77^ (£ > 0) of the order parameter El. If .-the order 
parameter is related to cluster multiplicity, like in the FIB processlafl, and if 
the most probable value of the order parameter is close enough to its average 
value, then (||) can be written in an equivalent form to the KNO scaling El . 
For large zm, the scaling function behaves as : 

$(z (1) )^exp(-z ( 1 1 / ) (1 - s) ) , (3) 

allowing for an alternative determination of anomalous dimension. For equi- 
librium systems at the critical point of second-order phase transition, g must 
be contained in between 1/2 and 1. 

We have to ask, what happens if the observable quantity is not an order 
parameter but the iV-dependent function of it like : m = N K — Nr), where 
K > g. For large N , \m\ is of order N K . Writing (||) with m instead of 77 and 
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Figure 1: The multiplicity distribution in FIB model for a = —1/3 with the Gaussian 
inactivation (ph (/3 = 0, c = 1, a = 1) in the scaling variables 8 = 0.76 are plotted for 
systems of different sizes : TV = 2 10 (circles), 2 12 (squares), 2 14 (crosses). 10 7 events have 
been calculated. Details of the calculation can be found in Ref. 11. 

taking into account : P[rj\drj = P[m]dm , one finds the generalized law : 

< |m| > 5 P[m] = $(z (a) ) = $ (^rjjj^j) , S= 9 -<1 (4) 

, which will be called in the following the 5 - scaling. The scaling function 
^(z^s)) depends only on the scaled variable : Zfg\ = (m — m*)/(< |m| > s ) . 
The scaling function ^(zrg)) i n (@) nas identical form as $(2/1)) except for the 
inversion of abscissa axe. In particular, its tail for large zi$) has the same form 

$(z (5) )~exp(-z ( 1 / (1 - 9) ) . (5) 

5 = 1/2 is a particular case when $>Qzrg\) is nearly Gaussian0 and < m >^ N. 
This limit, called the second ^scaling, has been found outside of the transition 
line in the shattering phase Ex More about the 5 - scaling interested reader 
will find in Ref. 9. 

The dynamical realization of S - scaling has been found in FIB model with 
finite-scale dissipation effects. For an appropriate choice of the fragmentation 



3 



kernel function, FIB is exactly equivalent to PQCD with the dissipation in 
jet fragmentation. In Fig. 1 we show the cluster multiplicity distributions 
calculated in FIB model for different initial sizes and plotted in the scaling 
variable z^ 76 y The calculation has been done for homogeneous fragmentation 
function : Fj^-j ~ [j{N — j)] a and the dissipation function which at small 
scales is approximated by the Gaussian inactivation rate function : 

I k = c^exp[-[(fe- l) 2 /N 2 ]/(2a 2 )]. (6) 

The parameters of fragmentation and inactivation kernel functions are : a — 
— 1/3, j3 = 0, c = 1, a = 1. Of course, the inactivation becomes scale-invariant 
if <T — ► oo and for this choice of a and j3, FIB process is in the shattering 
phase where the second - scaling holds. On the other hand, if a — ► 0, then 
Ik — > and the fragmentation process becomes independent of the cluster size 
k. This limiting situation happens at the transition line where the first-scaling 
holds. Thus, 5 - scaling in FIB appears rather as a cross-over phenomenon 
between first- and second- scaling domains than the phenomenon of gradually 
increasing deviation between the order parameter in the studied process and 
the relevant observable. Below, we shall address this question rigorously in the 
Landau-Ginzburg (L-G) theory of phase transitions. 

3 6 scaling in the Landau-Ginzburg theory 

Let us consider the homogeneous L-G free energy density as the simplest exam- 
ple of second-order phase transition. If e is the relative distance to the critical 
point, then the free-energy density is : f(rf) = erj 2 + brf + • • • , with b a positive 
constant. Under this form, the most probable value of ij is implicitely set to 
in the disordered phase. For finite size systems, it is often more convenient 
to work with the extensive order parameter : m = Nrj. Keeping the first two 
terms in f(rj), one may show that the probability that the system will be in a 
state to for a given e is : 

where is defined by the normalization of P[m]. Without loss of generality, 
we consider now the case when rj is contained between and 1, i.e., to > 0. 
We know that at the transition point (e = 0), the function P[fh] verifies the 
first-scalingDD , while below the critical point (e < 0), P[rh] verifies the second- 
scaling£3. These results are valid in the thermodynamic limit. How the finite 
system moves on between these two scalings when the control parameter e 
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moves away from 0? To answer this question, let us first write down the 
average value of m : 

< m >= iV 3/4 ^^ ^ '- . (8) 

f vr 1 ! 2 exp(— eN 1 / 2 u — bu 2 )du 

When N — > oo, then the term in brakets is only a function of eiV 1 / 2 . This is 
the finite-size scaling of this system near its critical point. It can be rewritten 
formally as : 

< m >= N 3 / 4 ip{eN^ 2 ) (9) 

At the transition point, ip has a constant value ip(0), an d < rh > behaves 
as a power of N with the mean-field exponent 3/4. At the critical point, 
this exponent is identical to the anomalous dimension defined in (Q). In the 
following discussion, it will be denoted by g both for systems at the critical 
state and outside of it. In the whole ordered phase, we have g = 1. Eq. (||) 
shows also that the cross-over between critical and non-critical phases must 
be discussed according to the values of eN 1 / 2 . Moreover, since we have to 
recover the non-critical case : < to >oc TV, for finite negative e therefore : 
ip(x) <~ (—a:) 1 / 2 , when x is large and negative. Let us note first the difference 
between the average value of to, as given by Q), and its most probable value 
(for e < 0) : rh* = ±N 3 l i [(-eN 1 ' 2 ) /2b] 1 / 2 . For e < 0, one can approximate 
P[rh] by a Gaussian function : 

N 3 / 2 l~4b / e n\ 

-P[rh] - J — exp(^-(TO-m*) 2 ) . (10) 



< TO > 



Suppose now that the order of magnitude of e is such that e ~ 1/N a , with 
some positive exponent a, smaller than 1/2 to ensure that eN 1 / 2 is a large 
number for large N . Then, < rh > and to* are both of similar order of 



magnitude. This means that the prefactor (TV 3 / 2 / < to >) in (10) is of order 
: < to >( 1+a )/( 2 ~ a ). In a similar way, the quantity (N/e) appearing in ( |ic| ) is 
of order : < to > 2 (!+ a )/( 2 - a ). Therefore, one can now write eq. ( |l0|) in the (5 
- scaling form (^). Obviously, we recover the known particular cases : S = 1 
for a — 1/2, and 6=1/2 for a = 0. Here, however, the ensealing appears as a 
finite-size effect and the tail of scaling function for large arguments has always 
Gaussian form, in contrast to the 5 - scaling due to N - dependent change of 
order parameter which was discussed in the preceding section. 

We have seen above that <5-scaling appears either as a result of incom- 
patibility between the observable and the order parameter , or as a cross-over 
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Figure 2: The UA5 multiplicity distributions in pp reactions □ for three different energies : 
200 GeV (circles), 546 GeV (squares) and 900 GeV (diamonds), are plotted in the S - scaling 
variables for S = 0.9. 

between first- and second-scaling. Both iV-dependent phenomena may actually 
be related one to another. In general, it is not so easy to conclude which of 
the two cases prevails without an additional information, and this point needs 
some comments here. Up to now, we did not need to know explicitely the 
total size N of the system, because the <5-scaling is not expressed explicitely in 
terms of N. In the practical applications, we have to know only that the data 
corresponds to a constant value of N. Nevertheless, we could have in principle 
access to the N - dependence of average quantities such as < m >. It should 
be : 

<m>~N 9,s , (11) 
for the change of variable as in Sect. || and : 

N 2g/(S+1) ; (12) 

for the finite-size cross-over effect described above in the mean-field approxi- 
mation. In the latter case, g = 3/4. If one knows from the experimental data 
the value of scaling parameter 5 and the A^-dependence of < m > , then one 
may find the value of the anomalous dimension if the cause of 5 - scaling can 
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Figure 3: The logarithm of scaled multiplicity distributions at UA5 dataf^l and three different 
energies : 200 GeV (circles), 546 GeV (squares) and 900 GeV (diamonds), plotted versus 

Z (0.9) = (( m " m *)/ < m >0 ' 9 ) 2 fOT Z (0-9) > °- 

be decided by analyzing the tail of the scaling distribution. Below, we shall 
illustrate this on the example of UA5 data. 

4 Application and conclusions 

We shall apply now the concepts developed above to the pp data of UA5 
Collaborations which show significant deviations with respect to the KNO 
(first-) scaling. If one plots data for different y/s in KNO variables, one notice 
that the maxima of < m > P[m] function are shifted , the height of maximum 
increases with y/s and the width decreases. Distributions at y/s — 200, 546 
and 900 GeV are plotted in Fig. 2 in the 6 - scaling variables for 6 = 0.9. 
Even if uncertainty on the average multiplicity is probably large, the scaling 
looks well. As noticed before, we can find argument for the cause of 6 - scaling 
by analyzing the tail of scaling function. This is shown in Fig. 3 where the 
logarithm of scaling distribution is plotted versus z^ Q g ^ . We see that the three 
curves tend to be linear in this plot, proving that tails of are essentially 

Gaussian. This provides a rather convincing argument that the fragmenting 
system is predominantly in the ordered phase (5 < 1, and Gaussian tail), close 
to a critical point {5 ~ 1). Moreover, a plot (unfortunately, for only three 



7 



values of y/s) of < m > vs y/s shows that : < m >~ (y/s) ' 35 . Using relation 
(|l2|), allows then to write : 

NO ~ (Vi) ' 33 , (13) 

what should be the proper scaling of the order parameter with the size of the 
system. Of course, one should be aware of the preliminary character of this 
extracted value of g. First of all, data in full space are not numerous. Secondly, 
the analysis is global in the sense that it concerns different classes of events 
mixed together. Nevertheless, the S - scaling analysis discussed in this work 
may become the powerful tool, allowing for an intelligent analysis of multiplic- 
ity data at different energies, which could give access to the determination of 
the anomalous dimension in the particle production at ultrarelativistic ener- 
gies. 
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